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Abstract. We compute the p-mode oscillation frequencies and frequency 
splittings that arise in a two-dimensional model of the Sun that contains toroidal 
magnetic fields in its interior. 



1. Introduction and Basic Equations 

Rotation and magnetic fields in the solar interior result in asphericity. This 
asphericity results in solar p-mode oscillation frequency splittings. From our 2D 
code we are able to obtain 2D solar models. In order to do helioseismic tests 
of the 2D effects, we need to calculate the p-mode frequency from the 2D solar 
models. The question is how to do so. 

No doubt, the starting points are the MHD equations. Nevertheless, it 
would be convenient to use the same format as in §19 of Unno et al. (1989). 
We reformulate them by defining the total pressure by adding the isotropic 
magnetic pressure component to the gas pressure, as we did in our ID and 2D 
solar variability model calculations (Lydon and Sofia 1995; Li et al. 2001, 2002, 
2003, 2006). The zero-order approximation of the MHD equations is equivalent 
to our stellar structure equations, which give rise to the equilibrium model of 
the Sun. We can obtain ID or 2D solar model by using our ID and 2D codes. 

The first-order approximation of the MHD equations yields the pulsation 
equations. The following adiabatic oscillation equations hold well for any equi- 
librium state when the magnetic fields are present, regardless its symmetry: 

v-e = -^(y + e-vp), 

a'^i + V0' + -Vp' + - V(/) = n', (1) 

where 

W = -^{(B- V)B' + (B'-V)B}, 

B' = V X X B). 

The symbols without primes represent the equilibrium state, which can either be 
spherically-symmetric, as for the ID stellar models, or aspherical, as for the 2D 
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stellar models that are to be extended in this paper. The symbols with primes 
are the Euler perturbations. It is well-known that represents the Lagrangian 
displacement vector. 

The adiabatic oscillation equations contain four equilibrium state variables. 
They are pressure (P), density (/)), gravitational potential (<I>) and the adiabatic 
index (Fi). All of them are functions of radial and colatitudinal coordinates (r, 9) 
in the azimuthal symmetry. We want to decompose them into two components 
in order to extend the oscillation equations into two dimensions like this: 

P{r,e) = P{r)+pi{r,e). (2) 

Here we assume that P{r) is the mean value of pressure over 9. As a result, 
pi{r,9) is the pressure asphericity. Similar for the others. This decomposition 
will automatically reduce to its ID counterpart when the equilibrium structure 
is spherically-symmetric. 

Substituting Eq. ([2|) into Eqs. ([T]) we obtain two-dimensional adiabatic os- 
cillation equations: 

1 a , 2. N 1 dlnP^ 

(r^^r) + + 



dr Ti dr 
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+ Vi) ^' - \-kGp - Ai^ = ^^GpC. 



The left hand sides are the same as their Id counterparts. The right hand sides 
represent the 2d effects, in which 
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In order to solve these equations we treat the two-dimensional correction 
terms like this: (1) to expand the oscillation variables in the 2D corrections 
in terms of spherical harmonics; (2) to remove the angular dependence of the 
2D corrections this way: (a) Multiplying them by the complex conjugate of the 
spherical harmonics; (b) Integrating them over the solid angle; (c) Summing 
them over all degree 1; (d) Summing them over all azimuthal order m. 
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2. Two-dimensional solar models and frequency splittings 



The torus field is defined as follows: -B^ = BqvWpo in the torus tube [{R — 
2a) sin O^ — r sin 6Y + [r cos 9 — {R— 2a) cos ^o] ^ = , but = outside the tube, 
where a is the tube radius, c= {R — 2a) sin is the distance of the tube center 
to the pole axis, is the colatitude of the tube center, R is the solar radius, po 
is the density at the solar surface, Bq is the magnetic strength parameter. 
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Figure 1. Dependence of density (left) and sound speed (right) asphericities 
on radius and colatitude in Model 1 

We want to investigate three 2D solar models: Model 1. a = 70 Mm, 
6*0 = 90°, Bo = 2400, B^rnax = 2.19 MG, resolution 2576 x 44; Model 2. The 
same as Model 1 except for resolution 2576 x 90; Model 3. The same as Model 
2 except for resolution 2576 x 180. The ID standard solar model (SSM) is also 
needed to compute the central frequency differences between 2D and ID models. 
Fig. [1] shows the density (left) and sound speed (right) asphericities in Model 1. 




1000 2000 3000 4000 5000 
u (mHz) 



X 

3. 









1 s e s 1 00 






101 S ( S 200 






201 S( S300 



































0.0 0.2 0.4 0.6 0.8 
Turning Point (r,/Ro) 



1.0 



Figure 2. Relative central frequency difference (left) of 2D model 1 with re- 
spect to SSM. Dependence of frequency splittings (right) on the inner turning 
points for Model 1 
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The central frequency is defined as u^i = i^nim=o, and the frequency sphtting 
is defined as Dnim = ^nim=o — ^nim- It is interesting to investigate how they 
change with the inner turning points. Fig. [2] shows the relative central frequency 
difference (left) of 2D model 1 with respect to SSM and the dependence of 
Frequency splittings on the inner turning points for Model 1 (right). Obviously, 
the splittings for whose turning points locate within the tube peak. This is easy 
to understand: their waves spend most time near their turning points since their 
group velocity vanishes at the turning points. 

Models 2 and 3 are designed to check how high the angular resolution is 
enough to minimize the numerical errors in the frequency splitting calculations. 
Fig.Oshows the splitting differences between Models 1 and 2 (left) and Models 2 
and 3 (right). Model 3 doubles the angular resolution of Model 2, which doubles 
that of Model 1. The fact that the splitting differences between Models 2 and 3 
decrease in comparison with those between Models 1 and 2 demonstrates that 
the angular resolution for Model 3 is not high enough yet. The main cause is 
that the step function gives rise to discontinuity. 




Turning Point (ri/Rg) Turning Point (rt/Rg) 

Figure 3. Dependence of splitting differences between Models 2 and 3 (left) 
and 3 and 4 (right) on the inner turning points. The shadow region is the 
location of the torus tube. 

We come to the following tentative conclusion: The frequency splittings 
change with the turning point this way: (a) the biggest when it is in the torus 
tube, (b) the smallest when it is above the tube, (c) in between when it is below 
the tube. We need to use a smooth field in order to obtain splittings of the 
required precision. 
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